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Nonselective quantum measurements, i.e., measurements without reading the results, are often 
considered as a resource for manipulating quantum systems. In this work, we investigate optimal 
acceleration of the Landau-Zener (LZ) transitions by non-selective quantum measurements. We use 
the measurements of a population of a diabatic state of the LZ system at certain time instants as 
control and find the optimal time instants which maximize the LZ transition. We find surprising non¬ 
monotonic behavior of the maximal transition probability with increase of the coupling parameter 
when the number of measurements is large. This transition probability gives an optimal approxi¬ 
mation to the fundamental quantum Zeno effect (which corresponds to continuous measurements) 
by a fixed number of discrete measurements. The difficulty for the analysis is that the transition 
probability as a function of time instants has a huge number of local maxima. We resolve this 
problem both analytically by asymptotic analysis and numerically by the development of efficient 
algorithms mainly based on the dynamic programming. The proposed numerical methods can be 
applied, besides this problem, to a wide class of measurement-based optimal control problems. 
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I. INTRODUCTION 

The control of atomic and molecular systems with 
quantum dynamics is an important branch of modern sci¬ 
ence with multiple existing and prospective applications 
in physics, chemistry, and quantum technologies [THU]- 
Significant efforts are directed towards the development 
of efficient methods for finding optimal controls for quan¬ 
tum systems |15H20) . In this work we approach this prob¬ 
lem for optimal acceleration of transitions in the two-level 
Landau-Zener system. 

The Landau-Zener (LZ) system [2TH24] is a simple 
model of nonadiabatic transition caused by avoided en¬ 
ergy level crossing. This system has a very wide range of 
applications in physics, chemistry, and biochemistry; see 
review in [25]. Some of the applications include trans¬ 
fer of charge |50|, photosynthesis chemical reac¬ 

tions [551 HH] 7 and manipulations with qubits. The last 
application includes the utilization of the Landau-Zener- 
Stiickelberg interferometry which is based on 

multiple traverses of the avoided crossing region, func¬ 
tional variation of time-dependent parameters of the 
model with a single traverse of the avoided crossing re¬ 
gion [HI HH |3M], and coupling the LZ system to an 
external environment |37H40| . 

In this paper, we use another quantum control 
paradigm—measurement-based quantum control, which 
was proposed and developed in |1TH15] . Experimental 
demonstrations of measurement-only state manipulation 
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on a nuclear spin qubit in diamond by adaptive partial 
measurements is described in [44] . The use of measure¬ 
ments as control is now actively studied also in com¬ 
bination with feedback [JS] |35] . Machine learning was 
used to generate autonomous adaptive feedback schemes 
for quantum information m- The measurement-based 
quantum control scheme [41] utilizes non-selective quan¬ 
tum measurements, i.e., measurements without reading 
the results, to manipulate quantum systems. Such mea¬ 
surements do not increase knowledge about the quantum 
state of the system, but, in contrast to classical mechan¬ 
ics, they influence the system and can be regarded as 
“kicks” destroying quantum coherence. We consider the 
problem of maximization of the LZ transition probability 
from one adiabatic energy level to another using a fixed 
number of non-selective measurements of adiabatic level 
at certain time instants. The goal is to find optimal time 
instants to maximize the transition. 

If the number of measurements is infinite and measure¬ 
ments are performed infinitely frequently, then the tran¬ 
sition probability tends to one due to the fundamental 
quantum Zeno effect |48H50j . However, an infinite num¬ 
ber of measurement is never experimentally achievable. 
So, one may ask a natural question of how to find optimal 
instants of a given finite and fixed number of measure¬ 
ments to approximate the quantum Zeno effect. Approx- 
imaton to the quantum Zeno effect by finite number of 
measurements was considered in |5I1455| without formu¬ 
lating the problem of optimality. The problem of optimal 
approximation was stated and solved for two-level sys¬ 
tems in im |43| in the case when measured observables 
can be arbitrarily chosen. In the present work, for a par¬ 
ticular model (the LZ model), we consider and solve an- 
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alytically and numerically the more restrictive case when 
observables are fixed while the measurement’s time in¬ 
stants are optimized. The advantage of this scheme is in 
its relative simplicity for experimental realization, since 
instants of measurements of a fixed observable (popula¬ 
tion) are simpler to modify than measured observables. 

We find a surprising nonmonotonic dependence of the 
maximal transition probability on the coupling strength 
when the number of measurements is large. Recently a 
nonmonotonic dependence on the coupling strength and 
temperature was found for a different physical situation 
of a LZ system coupled to a harmonic-oscillator mode 

EH- 

The following text is organized as follows. Section [ll| is 
devoted to a discussion of the model and the LZ formula 
for the transition probability. In Section |III| we discuss 
a measurement-based control scheme and formulate the 
optimal control problem. We show that the target func¬ 
tion, which is the transition probability as a function of 
time instants, may have a huge number of local max¬ 
ima. This feature makes it hard to solve the problem by 
well-known general global optimization methods such as 
random search, simulated annealing, evolutionary algo¬ 
rithms, etc. 

In Section |IV[ we develop an efficient numerical algo¬ 
rithm based on dynamic programming to solve this prob¬ 
lem. The algorithm is still computationally costly and in 
order to unravel mathematical structures standing be¬ 
hind this problem and to develop simpler algorithms, we 
consider limiting cases of antiadiabatic (Section]^ and 
adiabatic (Section |VI[ ) regimes. A graphical presentation 
of the results is given in Figs. and 


II. LANDAU-ZENER MODEL 


The LZ Hamiltonian at time t is defined as 

iL(t) = AtJx - etcTz, (I) 

where 



are Pauli matrices, and e and A are positive constants. 
The time evolution of a pure state \'4’{t)) of the system is 
defined by the Schrodinger equation (in the units h = 1) 

Time-dependent eigenvalues (adiabatic energy levels) 
and eigenvectors of the Hamiltonian have the form 

= tVA 2 -1- (£t)2, 


Energy Levels 



aj{t) G C. If aj{t) are such that |<^j (t)) are unit vectors 
for all t, then |(/3o(i)) and |(/ 7 i(—t)) tend to | 0 ) (up to a 
phase term) as t ^ -too and to |I) as t —>■ —oo. Here, 

|0) = (l O)^, |1) = (O l)^ is the standard basis of C^. 

The dependence of energy levels on time is shown on 
Fig.g A transition between the energy levels is possible 
mostly in the time interval when they are close to each 
other (so called avoided crossing region). 

A pure quantum state is specified by a unit vector 

IV’(O) = a(t)|0) -b5(t)|l), 

where a{t) and b{t) are complex coefficients, and |a(f)p-|- 
| 6 (t)P = 1. Expressing these coefficients in the polar 
form as a{t) = cos( 0 / 2 )e*‘^L b{t) = sin(d/ 2 )e*‘^^, 0 < 0 < 
TT, 0 < ipi ^2 < Stt, allows one to represent the state of the 
system as a point on the unit two-dimensional sphere, 
called the Bloch sphere, with spherical coordinates 9 and 

ip = (p2- Ifi- 

The Schrodinger equation takes the form 

{ id = —eta + Ab 
ib = Aa + etb. 

By the time rescaling r = V^t (and redenoting the 
new time variable t again by t) this system is reduced to 

{ id = -^a + 

. t P) 

lb = v^a-b 

where 7 = jle. This system of equations defines a 
family of unitary operators C/(t,to)j fojf G K- Let us 
denote their matrix elements in the standard basis as 

/uoo(t,to) Woi(t,to)\ 

\uio(t,to) uii{t,to)J ’ 


where j = 0,1, the top sign in ± and corresponds 
to j = 0, and the bottom sign corresponds to j = 1, 


where 


Uki{t,to) = {k\U{t,to)\l). 
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Then Moo(i,io) = a(t) and Mio(t,io) = ^(i)i where a{t) 
and b{t) are the solutions of system ([^ with the ini¬ 
tial conditions a(to) = 1, b{to) = 0. Analogously, 
Moi(t, to) = a(t) and uii(t,to) = t)(t) if the initial con¬ 
ditions are a{to) = 0, b{to) = 1- 

The general solution of system ^ can be expressed as 
(detailed derivation can be found, e.g., in [57]) 

a{t) = 

(3a) 

b{t) = 

(3b) 

Here D,^{z) is the parabolic cylinder function defined for 
arbitrary complex order i/ and argument z, D*{z) = 
Djj{z) is a complex conjugate function, and ci and C 2 
are arbitrary constants. The parabolic cylinder functions 
are related to the Whittaker functions and the confluent 
hypergeometric functions; see [58]. 

We will use the following asymptotics of the functions 
and D_iy_i{±e^'^^'^t) as t —>■ -l-oo: 

(4a) 

-p (4b) 

= 0(t"^), (4c) 

/rt —7r'*)'/4 

Z4_.,_i(-e*-/4t) = r(r+z7) 

Let us note that both functions are bounded as t € K. 
Consider the initial conditions such that 

|a(—oo)| = 1, b{—oo) = 0. (5) 

Asymptotics Q imply that ci can be choosen as 
g-’r7/4g3i7r/8 corresponds to the phase of a{t) equal 

to e®* for large negative t) and C 2 = 0, so that 

(6a) 

h{t) = (6b) 

This gives the celebrated Landau-Zener formula: 

|a(+oo)p=e-^-^ (7) 

The quantity |a(-|-oo)p is the probability that the final 

diabatic state of the system will be the same as the initial 

diabatic state. This probability will be of our interest in 
the following. Since it corresponds to a jump from one 
adiabatic energy level to another (the state |0) corre¬ 
sponds to large negative energy for t —)■ —oo and to large 
positive energy for t —>■ —oo), we will refer to it as “tran¬ 
sition probability” (sometimes the term “LZ transition” 
is used in the opposite sense: for a transition between 
diabatic states |0) and |1) and, so, staying on the same 
adiabatic energy level). 

As we see, the transition probability is close to zero 
for large 7, or, equivalently, for small s with respect to 


A^. That is, for very slow (adiabatic) evolution of the 
Hamiltonian, the system stays on its initial adiabatic en¬ 
ergy level. This limiting case is called adiabatic. In the 
opposite limit of small 7 , called the anti-adiabatic limit, 
the LZ transition probability is close to one. 

III. MEASUREMENT-BASED QUANTUM 
CONTROL 

A general state of a quantum system is specified by a 
density operator p which is a self-adjoint positive unit- 
trace operator. Its unitary evolution is given by 

p{t) = U{t,to)p{to)U'^{t,to), 

where unitary operators U{t,tQ) are defined above. 

A von Neumann observable for a two-level system is 
specified by two linear operators Pq and Pi such that 
PjPk = PjSjk, where j. A: = 0,1 and 6jk is the Kronecker 
symbol, and Pq + Pi = 7, where I is the identity operator 
(operator P satisfying the property P^ = P is called 
projector). 

We will use only non-selective measurements, i.e., mea¬ 
surements without reading the results. Such measure¬ 
ments do not increase our knowledge of the system’s 
state, but they cause an instantaneous change of the state 
to 

p' = PopPo + PipPi, ( 8 ) 

where p is the density matrix just before the measure¬ 
ment. In particular, if we perform a measurement in the 
standard basis, i.e., Pq = | 0 ) ( 0 |, Pi = | 1 ) ( 1 |, then 

p' = |0)(0|p|0)(0| + |l)(l|p|l)(l|. (9) 

Alternatively, one can represent the change of the quan¬ 
tum state as follows. Every two-dimensional density ma¬ 
trix in the standard basis can be represented as 

_ I ^ 1 + Wz Wx- 'iwy\ 

^ ~ 2 [wx + iwy 1-Wz ) ' 

where w = {wx,Wy,Wz) € is called the Bloch vector, 

I to I < 1. Measurement in the standard basis changes the 
state to 

, _ 1 I'l + Wz 0 \ 

^ ~ 2\ 0 1-WzJ- 

So, it eliminates the off-diagonal elements of the density 
matrix, or, equivalently, projects the Bloch vector onto 
the vertical axis. Unitary evolution does not change the 
length of the Bloch vector. 

In the density matrix formalism, initial conditions (§ 
have the form p(—oo) = |0)(0|. The LZ transition prob¬ 
ability, i.e., the left-hand side of Q, is now expressed as 
( 0 |p(-koo)| 0 ). 

We consider the following optimization problem: given 
natural JV, find instants of measurements ti < t2 < ■ ■ ■ < 
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tN such that the transition probability (0|p(+(X))|0) is 
maximal for the initial condition p{—oo) = |0)(0|. The 
system evolves according to the unitary evolution be¬ 
tween measurements and with jumps according to for¬ 
mula Q at the instants of measurements. 

The transition probability can be expressed as 

N 

(0|p(-|-oo)|0) = ^ ^ jfc^fc+i)I ) 

( 10 ) 

where to = —oo, tAr+i = -l-oo, and jo = Jn+i = 0. 

We could use 0 to express the matrix elements 
Ujk(t, to) by the parabolic cylinder functions and to solve 
the maximization problem numerically. However, the 
matrix elements are high-oscillating functions with re¬ 
spect to both to and t. Moreover, the frequency of os¬ 
cillations tends to infinity as 7 —> 00. This causes a 
huge number of local maxima in the corresponding con¬ 
trol landscape. As examples, target functions for some 
cases are presented on Figs. and Our attempts to 
apply the known methods of global optimization such as 
random search, simulated annealing, evolutionary algo¬ 
rithms, etc., show that the problem is hard to solve by 
these methods. 

First let us establish an upper bound for the maximal 
transition probability. 

Theorem 1. For all values of the arguments we have 


Further relations between these two measurement- 
based optimal control problems are considered in 
Sec. IVIBI 

Fig. [2] suggests that, for all 7, the optimal instant of 
the measurement in the case of = 1 is t = 0. We prove 
this rigorously for the limiting cases of small and large 
7. The corresponding value of the transition probability 
equals to 

(0|p(+oo)|0) = i(l + e-2-^). (12) 

This formula can be derived from 

0(//2 

^.(0) = (13) 

where r(a:) is the gamma function, which immediately 
follows from the definition of the parabolic cylinder func¬ 
tion and the confluent hypergeometric function, see |58j . 

This gives an important concluson that a single op¬ 
timal measurement decreases the difference between the 
transition probability and unity twice. In particular, for 
large 7, the LZ transition probability exponentially tends 
to zero, while just one optimal measurement increases it 
to 1/2. 

IV. DYNAMIC PROGRAMMING ALGORITHM 


(0|p(+oo)|0) < ^ |l-t ^cos-^^^^ |. (11) 

Here Aip denotes the angle between the Bloch vectors 
that correspond to the states |0) and U^+oo, —00) |0). 


Here we solve the optimal control problem using the 
dynamic programming paradigm [59]. Denote by /n(p, t) 
the maximal value of (0|p(-|-oo)|0) provided that the time 
t state of the system is 

p(t)=p|0)(0| + (l-p)|l)(l| (14) 


Proof. Suppose we can arbitrarily choose not only time 
instants of measurements but also observables. Then the 
state of the system after the measurement is changed ac¬ 
cording to ([^, where an observable {Po,Pi} is chosen 
by us and may be different for different measurements. 
Consider maximization of (0|p(-|-oo)|0) under these con¬ 
ditions. This problem was solved in mills]; the maxi¬ 
mal transition probability is given by the right-hand side 
of (11). Our problem has an additional restriction that 
the observables are fixed. The optimal value of the tar¬ 
get function in the problem with an additional restriction 
cannot exceed the optimal value of the target function for 
the problem without this restriction. □ 


Remark. Formally, in the problem considered in [411 
143] ■ the instants of measurements ti,...,tjv are fixed. 
But, in the case of variable observables, this is not a 
restriction, because the measurement of any observable 
{Pi, P2} at an instant t is equivalent to the measurement 
of the observable {P{, P2} at the instant t', where 

p; = uit',t)p,u\t\t), j = i,2. 


and n measurements are allowed to perform after the 
time instant t. If n = 0, then no more measurements 
are allowed and we have only unitary evolution from t to 
-boo which gives 

foiP,t) =p|woo(+oo,t)p -b (1 -p)|uoi(-boo,t)p. 

For n > 1, we have the following recurrent relation: 

fn{p, t) = max fn-i{p |moo (i', AP + (1 - p)\uqi{ f ,t)\^,t'). 

(15) 

The target function is (0|/9(-boo)|0) = /iv(0, —00). Thus, 
the problem of multidimensional maximization is reduced 
to a sequence of problems of one-dimensional maximiza¬ 
tion. 

For further simplification, note that /„(p, t) = fnip'■, t) 
if p—1/2 and p' — 1/2 have the same sign for any n > 0, so 
that only whether p<l/2orp>l/2 matters. Indeed, 
the Bloch vector corresponding to a density operator of 
form ( [l4| has the length |2p — 1| and is directed upwards 
if p >1/2 and downwards if p < 1/2. The target func¬ 
tion is (0|p(-boo)|0), so, the Bloch vector corresponding 






5 




FIG. 2. Graphs of the transition probability for one measurement at the instant t (left column) and for three measurements 
(right column) at the instants —t, 0 , and t, t > 0, for 7 = 0.1 (top curves), 7 = 1 (bottom thick curves), and 7 = 5 (bottom 
thin curves). 



FIG. 3. Graphs of the transition probability for two measurements at the instants ti,t 2 (left column) and for five measurements 
(right column) at the instants —t 2 , —ti, 0, fi, t 2 , 0 < fi < t 2 , for 7 = 0.1 (top plots), 7 = 1 (middle plots), and 7 = 5 (bottom 
plots). The ranges of ti and t 2 for the middle and bottom graphs are the same as for the top graphs. 
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to the final state p(+oo) should be as close as possible 
to the Bloch vector (0,0,1), which corresponds to the 
state |0) (0|. The evolution of a Bloch vector under the 
action of both unitary evolution and non-selective mea¬ 
surements is uniform: if a Bloch vector w(t) evolves into 
w{t'), then the Bloch vector Xw{t) evolves into Xw{t'), 
0 < A < Hence, the optimal control depends only 

on the direction of the initial Bloch vector, but not on 
its length. 

However, the optimal value of the target function de¬ 
pends also on the length of the initial Bloch vector. If 
p = 1/2, then the quantum state is p = 1/2 (sometimes 
it is called the chaotic state) and, by induction, one can 
show that fn{l/2,t) = 1/2 for all n and t. As p in¬ 
creases (decreases) from 1/2 to 1 (0), the function fn{p, t) 
changes linearly from 1/2 to [/„(0,t)]. So, 

f( /(l-2p)[/„(0,t)-l/2] + l/2, p<l/2, 

^ |(2p-l)[/„(l,t)-l/2] + l/2, p>l/2. 


and it is sufficient to solve the maximization problem in 
(15) only for p = 0 and p = 1. For a fast and precise 


solution to this problem we use the so called very fast 
simulated annealing algorithm. 

Now we can solve our maximization problem as follows. 
Successively, we calculate functions fnip, i)i P = 0,1, n = 
1, 2,..., for t S [—50, 50] with the step 0.01. This is 
sufficient for 7 < 5 and N < 15, i.e., all optimal instants 
of measurements certainly lie within this interval. The 
results of the calculations are presented in Fig. [^ and 

This algorithm is still computationally costly. In the 
next two sections we develop simpler algorithms and re¬ 
cover mathematical structures standing behind this prob¬ 
lem. Toward this aim, we analyze the limiting cases of 
small and large 7 corresponding to antiadiabatic and adi¬ 
abatic regimes. 


V. ANTI-ADIABATIC REGIME 


The anti-adiabatic regime is the case of small 7. It 
turns out (see the derivation in Appendix]^ that, in the 
first-order approximation with respect to 7, we have 


\uoo{t,-oo)\'^ 
|uoo(+00, to)P 

(0|p(-foo)|0) 


1 - ^ - TTjV2^{e^^/*T*{t/V2)} - jf^2F2{-ity2) + 

1 + ^{2TT-fT{to/V2)[T*{t/V2) - T*ito/V2)] - 2F2{-it^ 12) - 


■7^0 


7r7 


n{27Tj:F{to/V2)[^^ - F*{to/V2)] +jtl2F2{-itl/2)} 


2F2{-itl/2)} + 0{j^), 

+ 0(7"), 


N 

1 - 27rjY,^Wk/V2)[F*itk/V2) - F*itk+i/V2)]} + 0 ( 7 "), 

k^O 


(16) 


where <0 = —00 , Iat+i = -l-oo, and 3? stands for the real 
part. 

As we can see, we should keep only the terms with 
ji = j2 = ... = jAT = 0 in sum (10). Indeed, as we can 


see from (A2) and (A3), all of the other terms have the 
order 7^ or higher. 

The value of the Fresnel integral J^(t/'\/2) can be rep¬ 
resented graphically by the Kornu spiral. Thus, from the 
geometric point of view, the problem in the first-order 
approximation is to find points ti,... ,t]\[ on the Kornu 
spiral which maximize function (16) (see Fig. [^. 

Consider the case of a single measurement. Then the 
transition probability has the form 


(0|p(-|-oo)|0) = 1 - Try - 2Trj\J^{ti/V2)\'^ + 0(7^). 


Its maximum is achieved at = 0 and its value in the 
first-order approximation is 1 — Try. Since the first-order 
approximation to the LZ formula gives the value 1 — 2Tr7, 
this finding coincides with a conclusion in the end of 
Sec. [hi] that a single measurement reduces the difference 
between one and the transition probability twice. 


Taking the limit —>■ 00 in (16) gives the value one 

for the transition probability, wl 
quantum Zeno effect. 


rich corresponds to the 


We cannot solve maximization problem (16) analyti¬ 


cally for any N. This function also has many local max¬ 
ima for direct numerical solution. Now we develop a nu¬ 
merical method for solving maximization problem (16) 
using the dynamic programming. 

Let us define the functions (not to be confused with 
the functions /„ in Sec. IV) 


foit)=^\Tit/V2) 


F*it/V2)- 


g — iTr/i 


/„(!)= mm ^ 


/tv = min 


3? [F{t/V2) [F*{t/^/2) - F*{t'/^/2)^^ } 

^3? [F^/^^/2F*{t'/^/i) + 1 } + /Tv-i(i') 


where n = l,...,V—1. In the dynamic programming, 
one needs to consequently solve the minimization prob- 
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FIG. 4. Maximal transition probability vs number of measurements N for different 7 . Left (from top to bottom): 7 = 
0.01, 0.05, 0.1, 0.2, 0.5, 0.75. Right (from bottom to top): 7 = 0.75, 1, 1.2, 1.5, 2, 5. Inset: a part of the same graph. Results 
for 7 < 0.5 were obtained using the algorithm for small 7 , see Sec. [V] and verified by a more precise and general, but more 
computationally costly method of Sec. |IV[ Differences in the maximal transition probabilities found by these two algorithms are 
less than 0.01. Results for 7 > 0.5 were obtained using the general method of Sec. EYi but the differential evolution algorithm 
combined with the analysis of Sec. |VI| also gives good results. The difference in the maximal transition probabilities is less than 
0.01 for 7 > 1 and less than or of the order of 0.01 for 0.5 < 7 < 1. Th e curve for 7 = 5, as well as for larger 7 not shown here, 
almost coincides with the theoretical upper bound curve given by ( 20 l. 



N 



FIG. 5. Maximal transition probability vs 7 for a fixed num¬ 
ber of measurements; from bottom to top = 0,1, 3, 5,10,15. 
Gurves for = 0 and N = 1 are drawn using the exact LZ 
formula 0 and formula ( |12[ ) . Other calculations are the same 
as described in the caption of Fig. 


FIG. 6 . /n vs number of measurements N. The target func¬ 
tion is (0|p(-l-oo)|0) = 1 — 27 r 7 /iv-I- 0 ( 7 ^) [see (17l]. The case 
of no measurements corresponds to /o = 1 and the first-order 
approximation 1 — 27r7 of the LZ formula (for small 7 , or 
anti-adiabatic regime). 


lems for fn(t), n = 1,... ,N. The approximate value of 
the transition probability is then 

(0|p(+oo)|0) = l-27r7/^ + 0(72). (17) 

To solve the problem for arbitrary N, we successively 
calculate fn, n = l,...,Af, for t € [—10,10] with the 
step 0.01. This interval is large enough for Af < 15 since 
optimal time instants certainly lie within this interval. 

To improve the accuracy, we perform the following pro¬ 
cedure. Denote the solution found by the dynamic pro¬ 
gramming algorithm by ti,... ,t^. Then use the point 
(ti,...,tjv) in the Af-dimensional real search space as 


the initial point to search for the nearest local maximum 
of (16). This is much easier than the search for the global 
maximum. Denote the new point as (t^,..., t^y)- 

These operations do not depend on 7. According 
to (161, approximation of the target function depends 
on 7, but the optimal time instants in the first-order 
approximation do not. So, it is enough to perform these 
operations only once. The results are provided in Tablejlj 
Figures and are drawn using the data of Table jl] 

Our calculations show that these operations give very 
accurate results for 7 < 0.1; the value of the target func¬ 
tion at the point {t'l,... ,t'^) differs from the maximal 





























TABLE I. Solution of the problem in the first-order approximation with respect to 7 


A 

/o 

Optimal time instants 

1 

0 

= 0 

2 

0.188 

ti = -3.31, t2 = 0.12 

3 

0.360 

ti = -3.33, t2 = 0,t3 = 3.33 

4 

0.461 

ti = -3.38, t2 = -0.24, ts = 0.24, t4 = 3.38 

5 

0.521 

ti = -3.41, t2 = -0.39, ts = 0,t4 = 0.39, ts = 3.41 

6 

0.563 

= -3.44, t2 = -0.50, ts = -0.17, t4 = 0.17, ts = 0.50, U = 3.44 

7 

0.594 

ti = -3.46, t2 = -0.59, ts = -0.30, t4 = 0,t5 = 0.30, U = 0.59, tr = 3.46 

8 

0.619 

ti = -3.48, t2 = -0.66, ts = -0.39, t4 = -0.13, ts = 0.13, U = 0.39, tr = 0.66, tg = 3.48 

9 

0.640 

ti = -3.50, t2 = -0.73, tg = -0.48, t4 = -0.24, ts = 0, tg = 0.24,17 = 0.48, Ig = 0.73,19 = 3.50 

10 

0.657 

li = -3.52,12 = -0.78,13 = -0.56,14 = -0.34, Ig = -0.11, Ig = 0.11, h = 0.34, Ig = 0.56,19 = 0.78, lio = 3.52 

11 

0.672 

li = -3.54,12 = -0.83,13 = -0.62,14 = -0.42,15 = -0.21, Ig = 0, It = 0.21, Ig = 0.42,19 = 0.62, lio = 0.83, 

111 = 3.54 

12 

0.685 

li = -3.55,12 = -0.88,13 = -0.68,14 = -0.49, Is = -0.29, Ig = 0.10, ly = 0.10, Ig = 0.29,19 = 0.49, lio = 0.68, 
til — 0.88, ti2 = 3.55 

13 

0.697 

ti = -3.55, = -0.88, = -0.69, i4 = -0.49, is = -0.30, te = -0.10, tr = 0.10, is = 0.28, ig = 0.48, iio = 0.67, 

ill = 0.88,ii2 = 3.54,ii3 = 7.08 

14 

0.709 

li = -7.08,12 = -3.55,13 = -0.87,14 = -0.68,15 = -0.48, Ig = -0.29, I7 = -0.10, Ig = 0.10,19 = 0.29, lio = 0.48, 
111 = 0.68,112 = 0.87, li3 = 3.55, li4 = 7.08 

15 

0.720 

li = -7.09,12 = -3.56,13 = -0.91,14 = -0.73, Is = -0.55, Ig = -0.36, I7 = -0.18, Ig = 0,19 = 0.18, lio = 0.36, 

111 = 0.55,112 = 0.73, li3 = 0.91, li4 = 3.56, li5 = 7.09 



For larger 7 we perform one more procedure. We use 
the final point as an initial point to search 

for the nearest local maximum of the initial maximization 
problem (10) without assuming smallness of 7 . This gives 
very accurate in the above sense results for 7 up to 0.5. 
The results of this algorithm are presented on Fig.[^ 


An interesting finding is that the obtained optimal in¬ 
stants are symmetric with respect to zero except of the 
cases N = 2 and N = 13. This “symmetry breaking” 
occurs when a new spire for the Kornu spiral becomes 
populated, see Fig. But due to the symmetry of the 
problem, if ti,...,tjv is a solution of the optimization 
problem, then —tjq,... ,—ti is another solution which 
gives the same optimal value to the transition probability. 


FIG. 7. Points on the Kornu spiral that correspond to the 
optimal instants of measurements for N = 12 (sticks) and 
= 13 (balls), in the first-order approximation for small 7, 
or antiadiabatic regime. Points corresponding to ±00 are also 
marked. Note that the optimal instants of measurements are 
symmetric with respect to zero for A = 12 (this is the case for 
the majority of As), but they are nonsymmetric for A = 13. 
This “symmetry breaking” occurs when a new spire of the 
Kornu spiral becomes populated (see a point on a new spire 
in the top right corner of the graph) 


VI. ADIABATIC REGIME 

A. Analytic solution 

The adiabatic regime corresponds to large 7 . It turns 
out (see Appendix [b] for details) that, in this case, the 
transition probability can be approximately expressed as 


(0|p(-|-oo)|0) = -{1 — cosai ]^(cosai cosaj-i 

value of the target function found by the more precise *=2 

algorithm of Sec. |lV]by less than 0.01. -hsina; sinai_i cosyAw;) cosay}. 


( 18 ) 
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where ak = arccot(tfc/ 2 ^/ 7 ), k = 1,... ,N, 

uj{t) = sgn(i) ln(|i| + vT+T^) + t\/l+ 

Awfe = 01 ( 4 / 2 ^ 7 ) - w(4-i/ 2V7)- 
Let us consider the factors 

cos^AuJi = cos[7(wt, - 

Since 7 is a large parameter, we can assign to these fac¬ 
tors arbitrary values from [—1,1] changing the time ti by 
an infinitesimal value of order 2tt/juj'{ ti). This can be 
understood also from Fig. [^(bottom right, 7 = 5, fV = 3) 
which shows rapid oscillations and top and bottom en¬ 
velope curves; in the limit of infinitely rapid oscillations 
one can assign to the transition probability any value be¬ 
tween the two envelope curves. Assigning the value -1-1 
to all of these cosines gives 

1 f ^ 1 

(0|/9(-|-oo)|0) = 2 ) l~cosai cos(ai_i — a;) cosajv ^ • 

Put TT - ai = (fio, ai-i - ai = ipi-i (I = 2,..., A), 
aN = Then 


( 0 |p(+oo)| 0 ) = 


1 


N 


1 + cos ipi 
1=0 


and the initial maximization problem is reduced to the 
problem 


N 




1=0 

under the condition 


max 


N 




1=0 

It can be shown [331 solution is ipi = tt/{N + 1) 

for all 1. Hence, ai = tt — Trl/{N -|- 1), ti = — cos[Trl/{N + 
1)], ti = — 2^/7 coi{'Kl/{N 1)), / = 1,..., A. Thus, the 
optimal instants of measurements lie in neighborhoods of 
the instants 


4 = -2y^cot 


Trk 

N+1 


(19) 


with radii of the order of Att/ y/juj'(tk)- The maximal 
value of the transition probability is 


( 0 |p(+oo)| 0 ) = ^ I 1 + fcos 


iV+1' 


( 20 ) 


Thus, a particular choice of the values of cos yAw; gives a 
feasible solution. But, in view of Theorem[with Aip = 
TT in the right-hand side of 0 since, up to a phase, 
(-1-00,—00) jo) = |1) for large 7], it is impossible to 
get a solution that exceeds the right-hand side of ( 20 ). 
Hence, solution (19) is optimal. 


Remark. This optimal solution is not unique for A > 2. 
In the case N = 2, 

( 0 |p(+oo)| 0 ) 

1 


= 2 d - a[tiT2 + 


1 - tI){1 - t|) cos7Aa;2]T2 


(Tfe = cosafc). If we, as above, put COS 7 AW 2 = -|-1, 
then the solution is ti = — cos( 7 r/ 3 ) = —1/2, T 2 = 
— cos(27r/3) = 1/2. However, if we put cos 7 Aa ;2 = —1, 
but change the sign of, for example, ti, then, obviously, 
the value of the target function will be the same. So, 
in this case, ti = T 2 = 1/2. However, these are differ¬ 
ent time instants with an infinitesimal distance from each 
other, such that cos 7 Ati ;2 = — 1 - 

So, the analytic solution gives only neighborhoods of 
the optimal time instants. Precise values can be deter¬ 
mined numerically by finding the local maximum that is 
next to the point (4) • ■ • given by (19). Obviously, 
this problem is much simpler than a search for a global 
maximum in the initial problem. 

However, the difficulty is that the size of the interval 
47 r/^/ 7 a;'(tfe) indefinitely increases as 4 approaches zero, 
which happens for large N. If this size is comparable to 
or even larger than the differences | 4 - — 4 -±i|j then this 
approximation does not work. In this case, we propose 
the following approaches. 

The first approach is to find optimal instants that 
are close to zero by another numerical global optimiza¬ 
tion method, e.g., random search or simulated annealing 
methods. This problem is low-dimensional and can be 
solved efficiently. For larger times, the approximation of 
large 7 can still be used. We developed such a general¬ 
ization based on the results of Sec. IVI AI 

The second approach is to use the differential evolution 
algorithm (a kind of genetic algorithm or, more generally, 
evolutionary algorithm) to find the maximum in the not 
infinitesimal vicinity of the point (ti,...,tjv) given by 


(19). Our calculation shows that this gives good results 
such that the differences between the maximal transition 
probabilities found by this method and by the more pre¬ 
cise method of Sec |IV| are less than 0.01 for 7 > 1 and 
less than or of the order of 0.01 for 0.5 < 7 < 1. The 
use of the differential evolution method is much less com¬ 
putationally costly than the algorithm of Sec. |IV[ Recall 
that if 7 < 0.5, then the simple algorithm developed in 
Sec. jV] for the case of small 7 works well. Hence, in prin¬ 
ciple, for the whole range of 7 , we can use rather fast 
algorithms of Sec. jV] and of this section, instead of the 
computationally costly algorithm of Sec. EYl 

Let us stress that for the use of the differential evo¬ 
lution algorithm, it is important to restrict the search 
region to the vicinity of (4, ■ • ■, tiv) given by (19). If we 
search in the whole space K^, then the differential evo¬ 
lution as well as other algorithms of global optimization 
with the same computational power and time give large 
error: they find transition probabilities that are less than 
optimal by 0.1 or more. So, the use of the analytic results 
of this section is important. 
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B. Relation to measurement-based optimal control or, equivalently, 
problem with variable obserables 

U''{tk,to)PoU{tk,to) = Pwk, k= 


In [5TJ 133] the following measurement-based control 
problem for a two-level system is considered. Let there 
be no unitary evolution in the system; the initial state 
of the system be pure with the corresponding Bloch vec¬ 
tor Wq. We want to perform N non-selective measure¬ 
ments of arbitrary observables to maximize the probabil¬ 
ity of transition to another pure state with Bloch vector 
Wtarget ■ The angle between Wq and lOtarget is denoted by 
Aip. Let us describe the solution of this problem. Denote 
Wi,..., wn- as the vectors obtained by rotation of Wq by 
the angles kA(p/{N + 1) in the plane formed by Wq and 
Wtarget- Then the optimal observables are the projectors 
onto the states with the Bloch vectors Wi,... ,wn. De¬ 
note them as P^i, ■ ■ ■, P-wk- The optimal value of the 
target function is 


1 

2 




N + lJ 


AT-1-1 ' 


This solution can be generalized to the case with uni¬ 
tary dynamics, which was also considered in mi- Let the 
system evolve according to the family of unitary opera¬ 
tors U{t,to)- We want to find N observables measured at 
fixed time instants tg < ti < ... < t]\[ < T to maximize 
the transition probability to a desired state ptarget at a 
final time T. This problem is equivalent to the above 
problem without unitary evolution and with target state 
C/t (T, to)ptargetU{T, to). If , . . . , are the optimal 
projectors in the evolution-free problem, then the opti¬ 
mal projectors for the problem with evolution are 


Uitk,to)ProkU\tk,to), fc=l,...,W 


In our case, to = —oo, T = -foo, the initial state is 
p(-oo) = | 0 ) ( 0 |, the target state is ptarget = | 0 ) ( 0 |, 
and I7'l'(-|-oo,-oo)ptargetC/(-l-oo,-oo) = |1)(1|. Thus, 
A(p = TT. Further, in our problem we have the hxed 
observable Pq = | 0 ) ( 0 |j but we can choose the time in¬ 
stants fi,..., tjv. Hence, we must try to choose instants 
of measurements such that 


U{tk,to)P.u,^U\tk,to) = Po, k = l,...,N, (21) 


where Wk are defined as above. Of course, this is impossi¬ 
ble in general. However, let us consider the trajectory of 
the Bloch vector corresponding to the state (t, -oo) | 0 ) 
for —oo<t< -boo and large 7 . The vector starts from 
the north pole as t = —00 and finishes in the south pole 
as f = -boo. As we see from (B2) and (B3|, its trajectory 
on the Bloch sphere is an arc, on which the Bloch vectors 
corresponding to the optimal projectors lie. This allows 
the target function to achieve the upper bound ( 20 ) in 
the limit of large 7 . 

Thus, in this case, we can speak about a kind of du¬ 
ality of the problem with hxed instants of measurements 
and variable observables and the problem with hxed ob¬ 
servables and variable instants of measurements. 

Fig. shows the dependence of the maximal transition 
probability on 7 for various N. We see that for a large 
number of measurements, the dependence is nonmono¬ 
tonic (also we see it in Fig|^. The initial decrease of 
the transition probability with increasing 7 is in corre¬ 
spondence with the LZ formula. But, for 7 greater than 
some value (approximately from 0.7 to 0.9 depending on 
N), the probability increases up to the limit (20). This 
is caused by rapid oscillations, which allow one to choose 
the time instants such that the “ideal” conditions ( 21 ) 
are approximately satished with high precision. 


C. Maximin problem 


As shown above, the transition probability oscillates 
with high frequency for large 7 . This is due to the fac¬ 
tors cos( 7 Aa;i) in (18). A precise choice of time instants 
in neighborhoods with radii of orders in/ytjuj'(^tk) allows 
one to assign definite values to these factors. However, 
these results can be used in practice only if the oper¬ 
ation time of the measurement device is smaller than 
in/(tk). Otherwise, if the actual instants of mea¬ 
surement can not be controlled with the necessary accu¬ 
racy, then the final transition probability can be far from 
optimal due to oscillations. This can be seen in Fig. 
(bottom right) or Fig. (bottom right). 

So, it makes it important to consider a maximin prob¬ 
lem, that is, the problem of maximization of the tran¬ 
sition probability with the worst instants of measure¬ 
ments within the specified neighborhoods. In other 
words, for given ti,. .. ,tn, the device “chooses” values of 
cos( 7 Aa;i) that minimize the transition probability. The 
problem can be formulated as follows: 


( 0 |p(-boo)| 0 ) = 


max min 
Acl)2 , ••• 


1 — cos ai ]^(cos ai cos Q!/_i 
1=2 


- sin ai sin ai-i cos 7 Aw;) cos 
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Let us analyze this problem. Since the terms 
cos( 7 Aa;;) enter the target function linearly, the minimiz¬ 
ing values are ±1. Hence, the problem can be rewritten 
as 


( 0 |p(+oo)| 0 ) = 

. 1 

max min - 

Q;i>...>Q;iv signs 2 


( 22 ) 


N 


1 — COS ai cos(a/_i zb ai) cos / ■ 


1^2 


Solving this problem is equivalent to solving 


N 


min max cos a 1 

signs 


q cos(q;/_i ± ai) cos aAT. 
1=2 


(23) 


Theorem 2. All solutions of maximin problem (22 1 are 
as follows: 

(1) ai = 7 r/ 2 , all other ai are arbitrary not exceeding 


ai; 

(2) ajv = 7 r/ 2 , all other ai are arbitrary not less than 


UN', 

(3) every ai is either 0, 7 r/ 2 , or tt, and = 7 r /2 at least 
for one I and the order is satisfied: ai > a 2 > ■ ■ ■ > aN- 
In all cases, the transition probability is (0|p(-|-oo)|0) = 
1 / 2 . 


Proof. Obviously, the transition probability for all three 
cases (l)-(3) is equal to 1/2. Let us prove that these 
and only these arguments give the optimal solution. It 
is sufficient to show that for other values of a/, we can 
choose the signs in (23) such that the right-hand side of 
(23) is positive. 

Let both ax and un belong to either [0,7r/2) or 
(7r/2, tt]. Then, cosai cosaN > 0. Let us assign all signs 
to All differences a/_i — ai belong to the interval 

(—7r/2,7r/2) and therefore cos(q;/_i — a/) > 0. Then the 
right-hand side of (231 is positive. 

Now let Q!i G ( 7 ( 72 , tt] and un G [0,7r/2). Then, 
cosai cosajv < 0. If the solution is not of form 3), then 
there exist ai G (0,7r/2) U (7r/2, tt). For definiteness, sup¬ 
pose that there exists at least one ai G (0,7r/2). 

We have ai > a 2 ... > a^- Let k denote a number 
such that ak G [7r/2,7r], au+i G (0,7r/2). Let us assign 
the signs for Z < fc and Z > Zc -I- 2 to ” and the sign for 
Z = fc -I- 1 to “-I-”. Then, 


k N 

P cos(ai_i - ai) > 0 , cos(ai_i - a;) > 0 . 

/=1 Z = fc+2 


But ak + ak+i € ( 7 r/ 2 , 37 r/ 2 ), so cos(afc + a^+i) < 0. 
Combining these gives 

k 

cosai JJcos(a/_i - a/)cos(afc +afc+i) 

N 

X Y\ cos(a/_i — a/) cosa/v > 0. 

□ 


Thus, if the operation time of a measurement de¬ 
vice is comparable or exceeds the period of oscillations 
Att/ y/juj'(tk), then the optimal solution is a single mea¬ 
surement at the instant t = 0 which corresponds to 
a = 7r/2; repeated measurements at this instant as well 
as measurements at t = ±oo, which correspond to a = 0 
and a = TT, do not affect the transition probability. In 
this case, the transition probability has the form 

(0|p(+oo)|0) = ^^, 

It has no oscillating terms which otherwise would require 
a precise choice of instants of measurements. 


VII. CONCLUSIONS AND DISCUSSION 


This work analyzes the measurement-assisted acceler¬ 
ation of transitions in the LZ system. Control by nonse- 
lective measurements is considered as an alternative tool 
to other methods of controlling the LZ system. The main 
numerical method exploited in this work is dynamic pro¬ 
gramming which seems to be natural for measurement- 
based optimal control problems. 

We have obtained analytic behavior for small and large 
values of 7 and developed various numerical algorithms 
for solving the problem of measurement-assisted accel¬ 
eration of the LZ transition. The combination of these 
methods allows one to solve the problem for the whole 
range of 7 . The results are presented in Figs. and 
The range [0, 5] for 7 on these figures covers values of 7 
in experimental works, e.g., in [31 [33 ISO]- The increase 
of 7 above the value 7 = 5 almost does not change the re¬ 
sults due to achievement of the upper bound (20). Hence 


the results are complete for practical applications. 

We discover a surprising effect of non-monotonic de¬ 
pendence of the maximal transition probability on 7 for 
a large number of measurements. Decrease of the tran¬ 
sition probability for small 7 is in agreement with the 
LZ formula, but for large 7 , oscillation effects become 
essential and can be exploited to increase the transition 
probability. 

The obtained transition probability in the limit of an 
infinite number of measurements tends to one that recov¬ 
ers the quantum Zeno effect. Surprisingly, as can be seen 
from Fig. [^ the convergence to this limit can be very 
slow for the LZ system, especially for values of 7 in the 
intermediate range 0 . 2 - 2 . 0 . 

The optimization techniques of this work can be ap¬ 
plied to the measurement-assisted control of the dy¬ 
namics of physical quantum systems if the description 
by Hamiltonian ([^ is a suitable approximation. The 
main experimental challenge may be in performing non- 
selective measurements fast enough so that the duration 
of a measurement is much less than the intervals between 
the optimal time instants. The inability to perform mea¬ 
surements with good enough time precision can produce 
lower values of the observed transition probability. The 
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impossibility of arbitrarily fast measurement was also ar¬ 
gued to be a bottleneck for the quantum Zeno effect due 
to the time-energy uncertainty relation m- The advan¬ 
tage of the methods of Sections |V] and |Vl] is in their 
robustness due to the final local search which can find 
optimal solutions for slightly distinct models, e.g., if mea¬ 
surement imperfections are small. 

The analysis of Section |VI C| also provides an answer 
for the adiabatic regime. In this regime, if time instants 
of measurements are not well defined and it significantly 
affects the desired transition probability, then the opti¬ 
mal solution will be to perform only one measurement 
at the avoided crossing time instant. This measurement 
will increase the transition probability from almost zero 
to 1/2, which agrees with the analysis of the quantum 
Zeno effect in the LZ system based on weak measure¬ 
ments [62] and the analysis of the LZ dynamics under 
strong dephasing [63l |64]. Generally, the formalism of 
continuous or weak measurements can be used if the du¬ 
ration of measurements is comparable to or larger than 
the intervals between the optimal time instants [65ll69j . 


If the description of the physical system requires small 
corrections to the Hamiltonian 0 , then the methods of 
Secs. [Vj and VI still can be applied due to the robustness 
of the final local search step. Strong deviations which in¬ 
volve more general and complicated Hamiltonians with, 
for example, nonlinear dependence on time, coupling to 
an environment, presence of noise, etc., or even with not 
exactly known Hamiltonians, are beyond the scope of 


this work. However, the general approach of Sec. IV 


based on dynamic programming seems to be natural for 
measurement-based optimal control due to the discrete 
nature of control and we expect that it can be applied 
even to such general situations. 
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Appendix A: Approximations of the unitary matrix 
elements for the anti-adiabatic regime 


The aim of this appendix is to obtain formulas (16). At 


first, let us obtain approximations for matrix elements of 
operators U(t, to) in the zeroth order approximation with 
respect to 7 . Substitution of 7 = 0 into ([^ gives con¬ 
stant populations |a(t)p and |6(t)P [we can use formula 
or just solve system (0 with 7 = 0 ]. 
In particular, under the initial condition |a(— oo)| = 1, we 
have |a(-|-oo)| = 1. This corresponds to the zeroth-order 
approximation to the LZ formula. 

To obtain the next approximation, we use the follow¬ 
ing formula for the derivative of the parabolic cylinder 
function HQ]: 




d"f 


7=0 






Here 


nt) = \ - 




is the Fresnel integral, 


F 

P-^Q 


(Xi , . . . , (2^ 

61 ,... ,&„ 


= E 

n=0 


have jSS] 




l + V2e^^/'^T*{t/V2) 


Then the matrix elements Ujo(t,to), j = 0,1, can be 
written as 

\ *7/2 


(ai)n • • • iap)n ^ ^ Mjo(/7o) = u^o(t7o)e*(‘ 

ibi)n---{bg)un\ 


(A2a) 


is the generalized hypergeometric function and (a)„ = 
a(a -|- 1) • • • (fl + ^ ~ 1) is a Pochhammer symbol [(a)o = 
Ij. Often we will omit the parameters of the generalized 
hypergeometric function. 

Since 

^ + 0(y3/2), 

we do not need the derivative of if we 

are interested only in the first-order approximation. We 


uoo{t,to) = 1 + TTjJ^{to/V2)[J^*{t/V2) - J^*{to/V2)] 


*7 1 1 + 2 ^ 


(A2b) 


uio{t,to) = -iy/Try[T{t/V2) -T{to/V2)] + 0{-f''^XA2c) 

Since equations ([^ transform one into another by com¬ 
plex conjugation and substitution t —>■ —t, we have 

Uj\{t^to) = q( t, to)- 
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So, we have calculated the first-order approximations of 
matrix elements of unitary evolution with respect to 7 . 

The separation of the term [(1 -|- -I- 

is motivated by unboundedness of the function 
2 ^ 2 ( 1 ,1; 3/2, 2; —it^/2) as t —>■ ± 00 , which we will see 
in a moment. Now we can explain it as follows. The 
term in asymptotics Q is equal to one by the ab¬ 
solute value, but its first approximation with respect to 
7 is 1 -|- f 7 ln|t|. This is an unbounded function, which 
approximates incorrectly for large t. This approx¬ 
imation is adequate only for bounded |t|. We need to 
consider t —>■ ±00 as well, and the LZ formula is for¬ 
mulated for this case, hence we must separate this term. 
Otherwise, as we will see, the first-order approximation 
of the unitary matrix elements will be invalid for large |t| 
and we will not be able to reproduce the first-order ap¬ 
proximation even for the case of no measurements (the 
LZ formula). 

Since the function is undefined for t = 0, we sepa¬ 
rate the term + which is equivalent to when 

t —± 00 . 

Now let us perform the limits to ~oo and t — > -l-oo. 
We use the following asymptotic expansion for the gen¬ 
eralized hypergeometric function m- 


2 F 2 



mm 

m) 


[K22{z) -I- L22{ — z)] 


here function 2 Tb is understood as just an asymptotic 
series. Symbol ~ denotes that the right-hand side is an 
asymptotic series of the left-hand side function. Further, 
L 22 {—z) is defined by formula (3) on p. 15 of [75] (the 
formula given in m is valid only for non-integer posi¬ 
tive a). An algorithm for obtaining this formula is given 
in IZSj. We do not write the whole asymptotic series here 
because of its cumbersomeness. For our purposes, we 
need only its first term: 

L22{-z) = -^[\nz-y{^)] + ... = —^[ln4z-HC]-H... 

\J 'KZ Z \J 'KZ 

(the next term has the order Xjz^y where ipix) is the 
the digamma function, i.e., logarithmic derivative of the 
gamma function, ^{1/2) = —C — 21n2, and C « 0.577 
is the Euler (or Euler-Mascheroni) constant. Thus, 

The function (t^/ 2 ) 2 F 2 (—it^/2) logarithmically di¬ 
verges when |t| indefinitely increases. But the divergent 
term is compensated by the term —| ln(l-|-t^)/(l-|-t§) in 
the expressions for Moo(t,to) and uii(<,to)- This was the 
motivation to separate the term [(1 -|- t^)/(l -I- 

Also we use 



as \z\ —>■ 00 , — 37 r /2 < argz < tt/2. Here 

K 22 {z) = z'^e^ 2 Fo{b — a,c— a\z~^), v = l + a — b — c. 


lim 

t—^dzoo 


Fit) = ± 


1 -I- i 
2 


gi'7r/4 


Series diverges whenever p > q + 2 and z ^ 0, but So, 


uooit, - 00 ) = 1 - ^ [1 + V2F^m*itjV2) 
uoo{+oo,to) = 1-1- TrjF{to/V2)[ 




(A3a) 




^ - F*ito/V2)] + ^ 2 F 2 i- 'f) + ^{ln[2(l + 4)] + C} - ^ + 0(7^),(A3b) 


iti 


iiwit, -00) = e—+ V 2 e-*"/ 4 _^(VV 2 )] + 0(7"/"), 

nio(+oo,to) = - V2e-^^^^Fito/V2)] + Oi^), 

ftoo(-l-oo, - 00 ) = 1 - 7 r 7 -I- 0 ( 7 ^), Mio(-|-oo, - 00 ) = e“*’"/^\/27r7 -|- 0 ( 7 ^^^). 


(A3c) 

(A3d) 

(A3e) 


Form this, we can immediately derive (16 1 . Also, we see 
that 


|uoo(+oo, -oo)|^ = 1 - 27r7 -|- 0 ( 7 ^), 

so that the first-order approximation of the LZ formula 
with respect to 7 is reproduced correctly. 

Remark. Since we need only absolute values of ma¬ 
trix elements of the unitary evolution operator, our trick 


with separation of the term [(1 -|- t^)/(l -|- is 

not necessary. Indeed, in the first-order approximation, 
the quantities include only the real part of 

2 F 2 , which does not contain divergent terms. How¬ 
ever, our trick is essential for problems where not only 
the transition probability but also the phase is important 
(for example, for measurement-assisted Landau-Zener- 
Stiikelberg interferometry) or for use of higher-order ap¬ 
proximations. 
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Appendix B: Approximation of the transition 
probability for the adiabatic regime 


mula for the parabolic cylinder functions miZi]: 


The aim of this appendix is the derivation of ( 18 ). In 


the case of large 7, we use the following asymptotic for- 


= expj^ - y(l - Iny) 
= expj^ + y(l - In 7 ) 




(1 + 0 ( 7 -^)), 


(Bla) 


- ^+i7w(t)| (^ 1 +)(l + 0(7 ^)), (Bib) 


where 

uj{t) = sgn(t) ln(|t| + \/l + + ti/l + 

Formally, these asymptotics are valid only for t ^ 0 . 
However, note that the limits t —)■ ±0 of the above asymp¬ 
totic formulas for coincide with each 

other and with the large 7 asymptotics of Dij{ 0 ). Indeed, 
in view of ( 13 ), 


^7(0) -- ^expj^ - y(l-ln7) 


(here a b means lim a/b = 1), which coincides 

7—>-oo 


with (Bla) for t = 0 . Similarly, the limits t —>■ ±0 of the 
above asymptotic formula for 

coincide with each other and with the large 7 asymp¬ 
totics of 

7^D_i_iT,(0) - ^exp|^-t^(l-ln7)-^|. 


So, formulas (Bl) can be used for all non-negative t. 
The substitution of these asymptotics in ([^ gives 




+C2e 


-i7w(t) 


' 1 , 1 - 


t 


( 1 - 


VT+¥J’ 

t 


(B 2 a) 


Vl + t^ 


+ 0(7 ) 


where new constants are 

f 7r7 17 Sin 

7 = 7 exp I — ± —(1 - ln7) -f — 

Here and in the following, the top sign in ± and =p cor¬ 
responds to j = 0 and the bottom sign corresponds to 


j = 1 . In the following, for simplicity, we will omit the 
terms 0(7“^) and solve the problem in the principal ap¬ 
proximation with respect to 7. For the initial conditions 
0(2777) = 1) ^(277 7) = 0) the coefficients are 


Cj = 




111± 


^0 


(B 3 ) 


2 \ J 

For the initial conditions 0(277^0) = 0 , 6(277^0) = Ij 


Cj = 


\ 


i 1 . 


^0 


•\/l +1 


(B 4 ) 


Let us introduce a new time variable 
t 


T{t) = 


71+7 




(t = ±1 correspond to t = ±00). From (B 2 )-(B 4 ), we 
derive 

7 ^ 7 ( 2771 , 27 ^ 7 )!^ = ^(l±TT-o) 


2 V (l-'ro)(l-'r^)cos7Awl, 

where j,l = 0, 1, r = t(1), tq = r(lo), Aw* = uJt - uJto, 
and the “plus” and “minus” signs in ± correspond to 
I = j and I = 1 — j, accordingly. 

Let 277 7 be the initial time instant (it is equal to —00 
for our problem, but, for a while, we will consider a more 
general case with an arbitrary Iq), with N measurements 
performed at the time instants 277 7 , ■ ■ • , 277 ^Af- Let 
us define functions (again, not to be confused with the 
functions /„ in Secs. llVIandlVl) 


fjitN) = |uo7+00,iAf)P, 

fj{tN-k,tN-k+l, ■ ■ ■ jtiv) = 

N 

E n l^ji+i j'i ^01 ; 

jiv-fc+i l—N—k 
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k = 1, 2,..., where jN-k = j, jjv+i = 0, and ijv+i = By induction, it is not hard to prove that 
+00 [cf. ®]. Obviously, 

■ ■ ■ , tw) = 

|■WOj(^Af-fc,tAf-fc-l)|^/o {tN-k, ■ ■ ■ ,tAf) 

+ \uij{tN-k,'tN-k-l)\'^fl{tN-k, ■ ■ ■ ,tN)- 


N-k +1 


= -<1±tn TiTi^i + - t^){1 - t^_^)cos^Auji 


l=N 


(0|p(+oo)|0) =^^hm^/o^( 2 v^to, 2 V 7 ti,..., 2 v^tAr) = - 1 - n 11 + ^^/ {1 - t^){1 - t^_^) cos-fAui 


N 

n[’ 


tn 


/=2 


where Jill at' i® understood as a product in reverse or¬ 
der, Ti = T(ti), Aloi = uj{ti) — w(t/_i), and the product 
n ^2 defined to be equal to one if iV = 1 (so, in this 


case we have 1 — tiTat = 1 — rf in the curved brackets). 

Since t; G [—1,1], we can write t; = coso;, ai G [0,7r]. 
This substitution gives exactly (18). 
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